Transport in Almost Integrable Models: Perturbed Heisenberg Chains 
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The heat conductivity ^(r) of integrable models, like the one-dimensional spin-1/2 nearest- 
neighbor Heisenberg model, is infinite even at finite temperatures as a consequence of the con- 
servation laws associated with integrability. Small perturbations lead to finite but large transport 
coefficients which we calculate perturbatively using exact diagonalization and moment expansions. 
We show that there are two different classes of perturbations. While an interchain coupling of 
strength J± leads to k(T) oc 1/J± as expected from simple golden-rule arguments, we obtain a 
much larger k(T) oc l/J'* for a weak next-nearest-neighbor interaction J'. This can be explained 
by a new approximate conservation law of the J-J' Heisenberg chain. 
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The thermodynamic properties of many experimental 
systems (like, e.g., KCuFg or MEM-[TCNQ]2) are well 
described by a one-dimensional (ID) nearest-neighbor 
spin-1/2 Heisenberg model 1]. In such systems, mea- 
surements of specific heat or susceptibilities are in quan- 
titative agreement with exact results derived from the 
Bethe ansatz. The situation is different when trans- 
port is considered. As is typical for an integrable sys- 
tem 0, the heat conductivity of the Heisenberg chain 
is infinite at any temperature while experimentally 
measured transport coefficients are finite. In real mate- 
rials the unavoidable presence of (small) perturbations 
like longer range spin-spin interactions, interchain cou- 
plings, disorder or spin-phonon interactions, which break 
the integrability, are expected to render the heat conduc- 
tivity finite. One can therefore expect that conductivi- 
ties at finite temperature T > are singular functions 
of terms which break the integrability. This has to be 
contrasted with the behavior of thermodynamic quanti- 
ties and most other correlation functions which — at least 
for finite temperatures — vary smoothly as a function of 
small perturbations (assuming that no phase transitions 
are induced). Obviously the general question arises of 
how transport can be calculated in "almost integrable 
models", i.e., how strongly is the transport affected by 
small couplings which break the integrability. 

This question is not only important for systems well de- 
scribed by integrable Heisenberg or Hubbard models but 
is also of relevance for a much broader class of quasi- ID 
materials. The reason is that effective low-energy theo- 
ries in ID are notoriously integrable. For example, an ar- 
bitrarily complicated two-leg spin-ladder is, at low ener- 
gies, well described by an integrable Sine-Gordon model 
as long as the energy gap is much smaller than mi- 
croscopic energy scales like J. The term "well-described" 
implies again that the integrable model can be used for an 
accurate description of thermodynamics. To understand 
transport, however, one has to study again the effects of 
small perturbation (suppressed by powers of Ae/J) on 
transport. 



A further reason for our investigations is the general 
theoretical question of how singular are integrable mod- 
els and how are they affected by perturbations. While 
the analog question is well studied in classical systems 
with a small number of degrees of freedom, c.f. the fa- 
mous Kol'mogorov-Arnol'd-Moser theorem I4|, not much 
is known in many-particle quantum systems. 

Heat transport in spin chains has been subject of in- 
tense experimental 0, IE ^'^d theoretical [3 S IS 
(ill research in the recent past. Numerical studies 
on small systems at high temperatures 0, 0, have 
shown that non-integrable models — in contrast to inte- 
grable ones — have finite transport coefficients. However, 
the regime of small perturbations, which is probably the 
most relevant experimentally, is not easily accesible by 
those methods due to the singular nature of conductivi- 
ties near the integrable point. (More results are available 
for classical systems, see, e.g., 

Ellis El)- Analytical 
approaches which calculate transport at low T based on 
the analysis of slow modes within the memory matrix for- 
malism 0, are only valid for systems not too close 
to an integrable point, as the corresponding slow modes 
have been neglected. One motivation of the present work 
was actually the question whether these approximations 
are valid. 

We will consider ID spin-1/2 models with the Hamil- 
tonian 

H = Ho + Hi+ Hi,^ (1) 
consisting of an integrable part 

i 

and a small perturbation Hi or Hi^± which breaks in- 
tegrability. We will consider two kinds of perturbations. 
First, we take into account that in reality next-nearest 
neighbor spins (nnn) are also weakly coupled. 

i?i = J' ^ (Sf S'f^2 + SfSf_^2 + ^' St 31^2) 1 (3) 
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with g' = J'/ J ^ 1. Alternatively, we consider the weak 
exchange coupling between parallel spin chains. Intro- 
ducing chain indices a, (3, we use 



(a/3) t 



(4) 



with — J\_l J ^ 1. 

The heat current operator Jq = ji is obtained in 
the usual way from the continuity equation dthi — ji — 
ji+i, where hi is given by H = y^ .- hi. In the following 
we use a symmetrized version |l7j : 



J 



{q/3) 

for A = A' = 1, to obtain for ji, 

^ J2 [2Sti • (S^ X St+,) + g\3SU ~ 4S," + SS^^,)- 

a 

(StixSf+i)+.9±E(S?-i-S?+i)-(SrxSf)]+0(5'2). 

(a/3) 

In the integrable case, H — Hq, where the heat current 
is conserved, [H, Jq] = 0, the conductivity k{uj) devel- 
ops a singular contribution at zero frequency, ReK(w) oc 
S{uj). This singular behavior suggests considering a per- 
turbation theory for the inverse of k{uj), rather than k{ll!) 
itself. We therefore focus our attention on the scattering 
rate r(w) [3 defined by 



X 



(5) 



where x — P{{Jq\ Jq))oj=o is the static susceptibility 
of the heat current. In the integrable case V{lo) = 0, 
which reproduces the singularity. Accordingly, T{uj) will 
be small for small g and we can therefore expand the real 
part of Eq. in r at least for finite w. 



Tcj^Re k{uj) « Re r(tj) + ^^^^ 



(6) 



where the left-hand side can be evaluated via the Kubo 
formula for small g. To leading order in g we can there- 
fore express Rer(ti)) = 'Y^'^=2 ^"^^(a-') in terms of a cor- 
relation function of dtJq — i[H, Jq] to get 

ff^rsM^-Re/ dte'^-+'°^\[dtJQ{t),dtJQ])o. (7) 
^ Jo 

Since in the unperturbed system the heat current is con- 
served, its time derivative is proportional to the pertur- 
bation: Jq = 0{g) and the correlation function Q can 
be evaluated with respect to the unperturbed Heisenberg 
model. 
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FIG. 1: Leading order contribution to the scattering rate T{uj) 
from exact diagonalization of a 20-site Heisenberg chain (thin 
solid line) with weak nnn coupling J' (at T = oo). To this 
order the scattering rate vanishes at u = , implying anoma- 
lous transport . The same result is obtained if r2(aj) is recon- 
structed from A'^ moments (A'^ = 10. ..26) using the maximum 
entropy method. The inset shows that r2(aj = 0) — > for 
TV ~> oo. 



Equation (Q is well known in the context of the Mori- 
Zwanzig memory matrix formalism Igj. While IS 
exact in the limit y — s- at any finite frequency, it is 
important to note that this is not the case for to = 
where the expansion ^ can become singular. However, 
in our context it is sufficient to know that Eq. (0 gives a 
rigorous lower bound to the heat conductivity k(uj = 0) 
for small g. This has been shown many years ago for 
systems which can be described by a Boltzmann equa- 
tion by Belitz ^2^. The generalization of this result to 
almost integrable models will be presented in a forth- 
coming paper Systematic improvements of Q can 
also be calculated within the memory matrix formalism 

i[ilim,j9.]. 

First, we consider the Heisenberg chain with a weak 
and isotropic (A' = 1) nnn coupling. Figure ^ shows 
the leading order contribution T2{to) to the scattering 
rate determined from an evaluation of Eq. ((JJ) for large T 
using exact diagonalization. As similar physical quanti- 
ties (at large T) have been reported to show surpris- 
ingly large finite size effects (not observed in our case) we 
have also reconstructed T2 from an analytic calculation of 
its first 26 moments, J^^ ^w"r2(w) = ([^t"" Vg, Jq]), 
using a high-temperature expansion for an infinite sys- 
tem. We have used various methods to obtain T2{oj ~ 0) 
from these moments including a continued fraction ex- 
pansion, the Nickel method |23l | and the maximum en- 
tropy method (23| (see Fig.QJ. Although the curves differ 
depending on which method is used for reconstruction, 
all methods consistently show that T2{uJ —> 0) vanishes. 
Our exact diagonalization results also show that this is 
not an artifact of the T — > cx) limit as the limit is smooth 
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FIG. 2: Third and fourth order contributions to the scattering 
rate for various system sizes (see Fig.Ql, the first nonvanishing 
contribution being of order <;*. Note that finite-size effects are 
smaU. 



(see, e.g., Fig.^. 

We would like to emphasize that the vanishing of the 
scattering rate r(0) to lowest order is very surprising 
both formally and physically. Formally, one would ex- 
pect that any "generic" correlation function of type 
has a finite a; = limit at any finite temperature. Phys- 
ically, golden-rule arguments suggest that the breaking 
of integrability leads to a decay rate of the heat current 
of order J'^ . In the following we will first investigate the 
role of higher order corrections and then the influence of 
other terms which break integrability. 

Corrections to F up to order J'^ are derived starting 
from Eq. ©, where our lowest order result, F2, is used to 
determine the term of order F^. The dtJq-dtJq correla- 
tion function is then evaluated to order J'^ and J'^ using 
the wave functions and energies obtained from the exact 
diagonalization of Hq. The results are shown in Fig. 12 
Since ReF(ix') has to be positive and F2(0) = 0, it is not 
surprising that F3(0) also vanishes. F4(0), however, is 
clearly finite. We therefore conclude that the heat con- 
ductivity in the limit J' — > 0, A = A' = 1, has the form 



J' 



0.054(1) 



T'^J'^f{T/J) T^J'-i 



for T ^ 00, (8) 



where / is an (unknown) function of T/J only, with 
f{x 00) « 18.5 estimated from our exact diagonal- 
ization results shown in Fig. |3 Together with the ana- 
lytical explanation given below this is the main result of 
our Letter. 

We start with the observation that the time derivative 
of the heat current is linear in g as [Hq + gHi, Jq] — 
0{g). How can the naive golden-rule argument which 
suggests a decay rate proportional to g^ fail? This can 
happen if the presence of slow modes modifies the long- 



time behavior of the dtJp correlation function as dis- 
cussed, e.g., in [1,0,0. We therefore try to construct 
a new slow mode of the perturbed system Ho+gHi start- 
ing from the conserved heat current Jg of the integrable 
model Hq. Hence, we seek a solution Ji to the equation 



[Ho+gHi,Jn + gJi] = 0{g''). 



(9) 



As [Hq, Jo] = 0, we have to construct a Ji with 

[HoJi] = -[Hi,Jo]. (10) 

Before constructing Ji, we investigate the consequences 
of its existence for the correlator ((T)). With Jq = Jg+gJi 
we find 

-zJq = [gHi,Jo] + [Ho,gJi] + 0{g^), (11) 
= g[Ho,Ji-Ji] + Oig^). (12) 

As a consequence, the leading order contribution F2(w) 
to the scattering rate — by partial integration — may be 
written as F2(a;) = uj'^A{uj), where A{uj) is the (Ji — Ji) 
self correlator in the unperturbed system. We therefore 
conclude that k{uj = 0) diverges at least as 1/g^ if Ji 
exists. This trick of studying "readjusted" approximate 
conservation laws may well be useful for many other sys- 
tems with slow modes. 

We turn our attention to relation (|10|l . To find a solu- 
tion Ji we make the most general ansatz for it. Ji is a 
translationally invariant operator of finite range consist- 
ing of a linear combination of products of spin operators. 
By inserting the ansatz into Eq. H10() we obtain a system 
of linear equations for the unknown coefficients. This 
overdetermined system of equations turns out to have a 
solution in the case of an isotropic (A' = 1) nnn pertur- 
bation of the Heisenberg model with 



Ji = -g'J^Y.^S, 



-i + S,+2)-(S, xS,+3). (13) 



The explicit construction of Ji proves the absence of a 
J'^ contribution to the scattering rate as discussed above. 
Note that it is not possible to construct a Ji such that 
the commutator in Eq. Q is of order g^ rather than g^. 

While H13() can easily be generalized to the case of an 
anisotropic XXZ chain with A 7^ 1, no solution for Ji 
exists in the case of an anisotropic nnn perturbation with 
A' ^ 1. We therefore expect (and confirm numerically) 
that in the limit of small J' and small but finite (A' — 1) 



jyr^ 



_ 0.21(2) JVt2 
J'2(l - A'Yh(T/J) J'2(l- A')2 



for T — !■ CX3, 
(14) 

where h is an (unknown) function of T / J only, the value 
of which we can determine from the results shown in 
Fig. 121 in the limit T ^ 00 . This figure also shows the T 
dependence of k for T > J where we use Eq. ^ and x is 
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FIG. 3: Calculated heat conductivity of the anisotropic 
frustrated chain as a function of temperature for various 
anisotropics A = A'. Inset: leading order contribution to 
the scattering rate for the isotropic case (thin solid line) as 
well as with weak anisotropies for T ^ oo. 




FIG. 4: Leading order (in J±/J) contribution to the scat- 
tering rate of weakly coupled spin chains for T oo. 
The finite value at = leads to a conductivity k ~ 
0.091(3)jV(^^iT^) per chain where Z is the number of 
nearest-neighbor chains. 

calculated to order g'^ using exact diagonalization. Large 
finite size corrections prohibit calculations for T J 
within exact diagonalization. 

In many experimental systems we expect that the 
leading term which breaks integrability arises from a 
weak coupling J± of chains, Eq. 10} (or spin-phonon 
interactions [8j) . For this perturbation, Eq. H13|l has 
no solution and k ~ i/Jj_ can be evaluated at high 
temperatures from Q using exact diagonalization, see 
Fig. 0] Our value for the ladder in the limit J± 0, 
K 0.18J^/(J? T^), seems to be consistent with results 
of Zotos [iflj obtained for finite J± using Lanczos diago- 
nalization. 

To summarize, we have analyzed the heat transport 
in spin chains near the integrable point. In the pres- 



ence of a small next- nearest neighbor coupling J', which 
breaks integrability, one can construct a new approxi- 
mate conservation law. As a result, the heat conductivity 
remains extremely high, k ^ 1/J'^. For other perturba- 
tions like a weak inter-chain coupling J± this construc- 
tion is not possible and k ^ ^/ J±- Thereby we have 
shown that transport in "almost integrable models" de- 
pends not only quantitatively, but also qualitatively on 
the precise way in which integrability is destroyed. It 
would be interesting to study experimentally systems in 
which the strength of J' and J± can be varied systemat- 
ically, e.g., by chemical substitutions or by pressure. 
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